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We consider the violation of the Leggett-Garg inequality in electronic Mach-Zehnder inteferome- 
ters. This set-up has two distinct advantages over earlier quantum-transport proposals: firstly, the 
required correlation functions can be obtained without time-resolved measurements. Secondly, the 
geometry of an interferometer allows one to construct the correlation functions from ideal negative 
measurements, which addresses the non-invasiveness requirement of the Leggett-Garg inequality. 
We discuss two concrete realisations of these ideas: the first in quantum Hall edge-channels, the 
second in a double quantum dot interferometer. 
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Bell inequalities set bounds on the nature of the cor- 
relations between spatial-separated entities within local 
hidden variable theories^. In contrast, Leggett-Garg in- 
equalities (LGIs) set bounds on the temporal correlations 
of a single system 3 -^, and are derived under the assump- 
tions of macroscopic realism (MR) and non-invasive mea- 
surability (NIM)£ . 

Bell and Leggett-Garg inequalities are related in that 
their assumptions both imply the existence of a classi- 
cal probability distribution that determines experimental 
outcomes. The probability amplitudes of quantum me- 
chanics allow for violation of these inequalities: with Bell, 
the violation is due to entanglement between the two sys- 
tems; with Leggett-Garg, the violation occurs due to the 
superposition of system states and their collapse under 
measurement. 

The simplest LGI, henceforth referred to as the LGI, 
reads 



K = C 2 i + C 32 - C 3 i < 1, 



(1) 



where C a p = (Q(t a )Q(tp)) is the correlation function of 
the dichotomous variable Q = ±1 at times t a and tp. 
Since the first experimental violation 6 - of this inequality 
with weak measurements of a superconducting qubit, the 
Leggett-Garg inequality has been experimentally probed 
in systems as diverse as photons^—, defects in diamonds 
centers^, nuclear magnetic resonance^, and phosphorus 
impurities in silicon^. Whilst the subjects of these stud- 
ies may not be macroscopic, the LGI performs a useful 
role for microscopic systems as an indicator that the de- 
vice is operating beyond classical probability laws. More- 
over, if one accepts that the alternative to classical proba- 
bilities is quantum mechanics, the LGI provides a decisive 
indicator of the "quantumness" on a system^ 3 .. 

In this paper, we are interested in the violation of the 
LGI in quantum transport, and in particular, in electron- 
interferometers. Although there has been much work 
on Bell inequalities in electron transport, e.g. Refs H^ - 
I22I the LGI has only relatively recently been consid- 
ered in this settin g 23 ' 24 . Specifically, the charge flow- 
ing through a confined nanostructure, e.g. double quan- 



tum dot (DQD), has been shown to violate an inequal- 
ity similar to Eq. (Q]) out of equilibrium^ 3 .. Further- 
more, the moment-generating function of charge trans- 
ferred through a device has also been shown to be subject 
to a set of LG-style inequalities, which are violated for 
various quantum dot models. The violation of LGIs in 
excitonic transport has also attracted recent interes t 25 ' 26 . 

There are several difficulties which make the investi- 
gation of the LGI in electronic transport challenging in 
practice. Ostensibly, the measurement of Eq. ([TJ requires 
time-resolved measurements where the time between suc- 
cessive measurements is smaller than the decoherence 
time of the system. For the double quantum dot of 
Ref. for example, this decoherence time is of the order 
of lna21, which makes the necessary time-resolved mea- 
surements very challenging (but, in principle, possible^). 
Furthermore, for the violation of Eq. dTJ to be a mean- 
ingful indicator of non-classical behaviour, it must be 
ensured that the measurements are non-invasive. This 
"clumsiness loophole"— that allows violations of Eq. ([1]) 
to be associated with invasiveness of measurement, along 
with possible circumventions, have been the subject of 
much discussion ^ 26 ' 30 . 

The transport set-ups we consider here are based on 
the electronic Mach-Zehnder Interferometer (MZI), and 
can overcome both of these problems. The basic idea is 
that an electron travelling through a MZI can take one of 
two paths, and this path index defines the variable Q = 
±1. Unidirectional passage of the electron through the 
system allows us to map the time indices of Eq. (JTJ onto 
positions within the interferometer. As we show below, 
this removes the need for time-resolved measurements. 

We consider two realisations of the MZI in which 
measurements of Q arc performed in two different 
ways. In the first, the MZI is formed from quantum 
Hall edge channels, a set-up which has been realised 
experimentally 3 ^— and also attracted a large degree of 
theoretical attention 3 ^— . By interrupting the edge 
channels at various points and diverting electron flow to 
current meters, we show that K can be obtained from 
measurements of mean currents alone. Furthermore, due 
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to the spatial separation of the Q = +1 and Q = — 1 
channels, our detectors interact with only one of the two 
Q-states at any given time. Thus, our scheme provides a 
natural way to implement ideal negative measurements, 
as advanced by Leggett and Garg as a way to satisfy the 
NIM criterion^. 

The second set-up we consider is a MZI with a quan- 
tum dot (QD) in each arm. This geometry is similar to 
several experiments^— that have investigated transport 
through Aharanov-Bohm rings with QDs in the arms. 
The difference here being that the dots are fed by two 
tunnel-coupled leads^i, rather than just one. The QDs 
are monitored by quantum point contacts, whose trans- 
mission is sensitive to the charge state of the QD^Sr— . In 
contrast with the first set-up, electrons are not diverted 
out of the MZI at any point, and the influence of the de- 
tectors occurs as a pure dephasing effect. The three cor- 
relation function in Eq. (TTJ) are obtained through a com- 
bination of mean currents, both through the MZI and the 
quantum point contacts, and zero-frequency noise mea- 
surements, which cross-correlate current fluctuations in 
the MZI and quantum point contacts. As in the previ- 
ous scheme, we construct an ideal negative measurement 
scheme with this set-up. 

Both of these techniques exploit a combination of su- 
perpositions of paths through an interferometer com- 
bined with a gathering of "which-way" information to 
violate the LGI. The first set-up is a particularly simple 
realisation of the LGI, and is by no means restricted to 
transport, but could be used e.g. with photons, atoms or 
molecules. 

The paper proceeds as follows. In Sec.UJwc describe the 
basics of testing the LGI in a MZI. Sec. |TT] describes how 
this may be translated into experiments with quantum 
Hall effect edge-channels. Finally, Sec. Mil considers the 
alternative DQD-QPC geometry studied here. 



I. MACH-ZEHNDER INTERFEROMETER 

We begin by describing an abstract version of our MZI 
scheme to outline the basic ideas. The MZI is a two- 
channel interferometer with two beam-splitters that di- 
vide the MZI into three zones which we label: 1, the 
input ports; 2, the arms of the interferometer; and 3, the 
output ports (see Fig. [TJ). We inject one electron at a 
time into the MZI and the path taken by the electron 
will be the degrcc-of- freedom under test with Q = +1 
when the electron is located in the upper channel of the 
MZI, and Q = — 1 the lower. Since the electron passes 
sequentially through the three zones, we can map a mea- 
surement of Q at time t a to a "which-way" measurement 
at any point in the region a of the interferometer. In 
particular, Qi and Q3 are measured at the input and 
output ports, and Q2 is measured by placing detectors in 
the arms of the interferometer at 2±, where 2 refers to the 
zone, and ± the upper or lower channel. In this section 
we assume that we have ideal single-electron detectors 



(a) 




p 

/ 3+ 




\j 2 




\ 3 " 
• 
















< 



FIG. 1. (Color online) The Mach-Zehnder interferome- 
ter with three different detector configurations for the non- 
invasive measurement of the LGI of Eq. fTJ). Electrons are 
injected into the 1+ port, (a) Complete MZI configuration 
with detectors only at the final outputs 3±. With this set-up 
we can measure the probabilities Pjy_(l) and construct C31. 
(b) An additional detector is inserted into the MZI + arm. 
With this configuration we can measure probabilities P£+(0) 
and P3±;2+(L ■)• ( c ) A detector in the ' — ' arm allows us to 
obtain P£L(0) and Pj± ; 2-(1, ')• Combining the results of (b) 
and (c) allows us to construct correlation functions C21 and 
C32. 



that "click" on detecting an electron, which is then re- 
moved from the system (i.e., the detectors act essentially 
as electronic analogues of photodetectors). More realis- 
tic measurements in terms of currents are discussed in 
section [TTJ 



A. Ideal negative measurements 

A detector placed in one of the arms interacts strongly 
with electrons in that path (they are completely removed 
from the MZI) and has no effect on electrons in the other. 
With a detector placed at Q = +1, say, then the ab- 
sence of a detector response (combined with MR and 
ideal detectors) allows us to infer the state of the sys- 
tem {Q = —1) without any disruption. This is exactly 
the form of detector required to perform an ideal negative 
measurement as envisioned in Ref. 0. 

To make the measurement scheme as simple as possi- 
ble, let us inject electrons into the 1+ port, such that 
the initial state is known^i. We do not need to mea- 
sure in zone 1 and there is no question about the NIM of 
Q\. The correlation function C21 and C31 boil down to 
measuring (Q2) and (Q3) respectively. 

Let us define P®± (n) as the probability that the detec- 
tor placed at position a± either detects (n = 1) or fails 
to detect (n = 0) the electron. Since no further measure- 
ments are made past point 3, it is irrelevant whether we 
measure non-invasively or not at point 3. Placing detec- 
tors at 3±, we measure the probabilities P£±(i), and the 
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C31 correlation function can simplify be expressed as 

C 31 =P 3 D + (1)-P 3 D _(1). (2) 

The set-up for this measurement is shown in Fig. [TJi. 

Since, in measuring C 2 i, no further measurements are 
made after region 2, it is also not necessary to measure 
C21 non-invasively. We can measure (Q2) (and thus C21) 
by running the experiment once with a detector in chan- 
nel 2+, and once in 2— (Fig. [T]d and c) and writing 

C 21 =P 2 D + (1)-P 2 D _(1). (3) 

It is perhaps instructive to discuss how to obtain this 
quantity using the ideal negative measurement technique 
and measure C21 in terms of the probabilities of absence 
of detector clicks. With the detector at 2+, we can equate 
the probability that no electron is detected, P£+(0), with 
the probability that the electron travels the path 2—. 
Swapping the detector to the other arm, we measure 
P^-iO) and infer the probability that the electron takes 
path 2+. Whence, we obtain the non-invasively mea- 
sured 

C 21 =P 2 D _(0)-P 2 D + (0). (4) 

Since P 2 D ±(0) = 1 - p £±0-)> E q- © and E q- © S ivc thc 
same result. 

We now consider C32 , where it is essential that we mea- 
sure Q2 non-invasively, since a subsequent measurement 
is performed. On the face of it, measuring C32 requires 
a correlation measurement between two detectors. This, 
however, is not the case, as we now show. 

Let us begin by placing one detector at 2+ and another 
one at 3+ (Fig. [T}d). We can then obtain the four prob- 
abilities, P^ 2 , (n, n'), that the detectors at 3+ and 2+ 
give the results n, n' = 0,1 respectively Of these, the one 
we are interested in is P® + 2 , (1,0), since this allows us 
to infer (non-invasively) the probability that the electron 
took path 2— to detector 3+. Moreover, we do not actu- 
ally need to actively detect at 2+, since, if the electron 
reaches the 3+ detector, it is clear that it has not en- 
tered channel 2+ (because the detector there would have 
removed the electron from the system)^. With all four 
probabilities, P^ q 23 /(l, ■), obtained in this non-invasive 
way, we can construct 

c 32 = - rfP&M 1 '-)* ( 5 ) 

where we have replaced the measurement value at posi- 
tion 2 with a dot to indicate that we do not actually have 
to measure there (the value is guaranteed to be zero). 

In this way we obtain all the required correlation 
functions, measured non-invasively where necessary. Al- 
though we have concentrated on the simplest case here, 
the above non-invasive techniques are extensible to the 
case where the input state is unknown and all C a p must 
be measured in a non-invasive way, or to more compli- 
cated LGIaii. 



B. Leggett-Garg Inequality 

The action of the MZI can be specified by two beam- 
splitter scattering matrices sx; X — A, B. With a aq the 
annihilation operator for an electron in channel aq, the 
beamsplitter input-output relations read 




Parameterizing the scattering matrices as 

/ cos(i^) sin(i^) e ^\ 
' X V-sin(i0 x )e- i ^x cos(i0 x ) J ' 1 ' 

we obtain the correlation functions 

C 21 =cos0 A ; (8) 
C31 = cos 6a cos Ob — sin 6a sin 6b cos 4>; (9) 
C 32 = cosd B , (10) 

such that the LG correlator reads 

K(6a, 6b, </>) = cos 9a + cos9b — cosOa cosOb 

+ sin6A sin 6b cos <f>, (11) 

with (f> = 2-((f>A ~ 4>b) being the phase difference accu- 
mulated between the two paths. This is a familiar ex- 
pression. If we identify 6a = and 6b = ^t 2 , then 
Eq. (|lll) is exactly that obtained for a qubit evolving un- 
der the Hamiltonian H = \£la x measured in the a z basis 
at times ti, t% = ii+Ti, and t 3 — t 2 + r 2 . Thc properties 
of Eq. (fTTj) arc discussed in Sec. UU 

II. QUANTUM HALL EDGE-CHANNELS 

Quantum Hall edge channels have been shown to al- 
low a direct translation of the MZI into electronic trans- 
port experiments^— and Fig. [2] shows a sketch of the 
quantum Hall geometry needed to realise our proposal. 
Each channel in the MZI is realised with a single edge- 
channel and the electronic beam-splitters are realised by 
quantum point contacts (QPCs). Backscattcring is sup- 
pressed between edge-channels such that transport is uni- 
directional. This set-up is the same as the MZIs of ex- 
periment except for the addition of extra contacts to the 
arms of the interferometer. These contacts are connected 
to the edge-channels via adjustable quantum point con- 
tacts, such that thc detectors can be coupled into and 
out of the MZI as required. This method of coup ling 
probes to the MZI arms has been realised in Ref. |38|. 
Port 1+ is raised to a voltage +V and electrons are in- 
jected into this channel. The output ports (detectors) 
are all grounded. When the correlation function C31 is 
being measured, the detectors at 2± are not required and 
are isolated from the MZI by closing their QPCs. (Fig. [2] 
shows detector 2— closed off in this way). To measure the 
remaining correlation functions, the detectors at 2± are, 
one then the other, connected into the MZI by opening 
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FIG. 2. (Color online) Quantum Hall edge-channel real- 
isation of the MZI set-up for measurement of Legget-Garg 
inequality. The MZI set-up is similar to that of Ref. [3l] but 
with two extra detectors (2±). These additional detectors can 
be isolated from the circuit by closing off the QPCs between 
them and the edge channel. The configuration shown has the 
detector at the 2+ position active such that transmission to 
beamsplitter B via channel 2+ is blocked, and the detector at 
2— is pinched off. This detector combination corresponds to 
that of Fig. [lp. 



up their respective QPCs. In Fig. [2 the detector at 2+ 
is connected into the circuit and fully prevents electrons 
in channel 2+ from reaching the outputs 3±. 



A. Current measurements 

Let (I aq ) be the mean stationary current flowing into 
output aq, given that when a = 3, the detectors at po- 
sitions 2± are closed off. Further, let (l3 q -2g'} be the 
current flowing at output 3g when the output at 2q' is 
open. Since, in the linear regime, the current operator 
for each output is I aq = GoVa^aag, with Go = e 2 /h 
the conduction quantum^, these mean currents are pro- 
portional to the probability that an electron travels in 
the corresponding channel. The correlation functions re- 
quired for the LGI can then be constructed, as with the 
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Division by the sum of detector currents removes pro- 
portionality factors and, if all detector are identical, also 
removes detector inefficiencies. Writing the scattering 
matrices as 



sx 



rx t' x 
tx r' x 



(14) 




FIG. 3. (Color online) (a) LG correlator K(0 A , 9 B , <f) of 
Eq. (|11[) as a function of the beamsplitter angle 9 = 9 a = Ob 
for three values of the phase (f> — 0, 7r/4, 7r/2. The shaded blue 
region indicated violation of the LG inequality (K > 1). The 
maximum violation, K max = |, occurs for = and, e.g., 
9 — tt/3. (b) The influence of dephasing. Shown is the LG 
correlator K of Eq. (|16p maximized over 9 a/ b as a function 
of the phase <f>. Results shown for values of the dephasing 
parameter A/tv = 0, |, 1, |, 2. Violations of the LG are only 
observed for cos0 > 0. 



we obtain the correlation functions 



C 2 i = \r A \ 2 -\t A \ 2 ; 
C 3 i = \r B rA + t' B t A \ 2 



\tBTA 



C 32 = \ta\ 2 {\tb\ 2 - \t B \ 2 } - \t A \ 2 {\t' B \ 2 - \r' B \ 2 } .(15) 

With scattering matrices as in the previous section, the 
LG parameter K obtained from current measurements is 
the same as Eq. (|11[) . This quantity is plotted in Fig.[3^i. 
A maximum violation of if max = § is obtained for pa- 
rameters 9 a = 9 B = 7r/3 and <f> = 0. 

The violation the LGI in this set-up arises because 
the measurements at 2± remove electrons from the in- 
terferometer arms, preventing interference between the 
two paths. The presence of this interference in C31 com- 
bined with its absence in C32 leads to the violation. 



B. Dephasing 

We can account for the effects of dephasing by allowing 
the phase <f> to fluctuate. We replace <f> — > 4> + Sip in 
Eq. (|11[) and integrate &p over a flat distribution in the 
range —A/2 < 8cf> < A/2. The resulting LG parameter 
with dephasing reads 

K dcph = cos 9 A + cos 9 B - cos 9 A cos 9 B 

+/(A) sin 9a sin 6 B cos 0, (16) 

with /(A) = 2A _1 sin(A/2) the function containing the 
dephasing effects^. If all angles are freely variable then 
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the maximum of this function is 

obtained for cos 0,4 = cos 0b = [1 + /(A)] -1 . Expanding 
for small A, we find if^ c a p x h (A) = § - ^A 2 . In the op- 
posite limit, where the dephasing is total, A — > 27r, we 
have /(A) — > and the maximised Leggett-Garg corre- 
lator reverts to the classical value, lim^^n ^2 = 1 as 
required. 

One interesting feature occurs if we assume that the 
phase <f> is fixed (e.g., we are not able to vary the magnetic 
field) and maximise over 6a/b (see Fig. [3b). Providing 
that coscj) > 0, the maximum value is 



K dcph tR J 

max(9 A/B )^ 



b,A) = /(A)cos0 + 



1 



1 + /(A) cost 



■(18) 



found by setting cos 9 a = cos 0s = [1 + /(A)cos(/>] 1 . 
If, however, coS(p < 0, the maximum value is just the 
classical value, KmaL(e A B ) = ^' f° un d by setting cos 6a = 
cos 9b = 1 • This reversion to the classical value occurs 
when the scalar product between the axis of the rotation 
of beamsplitter B and that of beamsplitter A becomes 
negative. 



C. Multi-channel case 

The above scheme is easily modified to include multi- 
ple channels. We take the same geometry as before but 
assume that each lead supports M channels. The M 
channels of the upper lead are all associated with qubit 
state Q = +1; the M channels in the lower lead, with 
state Q = — 1. The scattering matrices of Eq. (|14p are 
thus generalised to 2Mx2M matrices with MxM blocks, 
Tx-, tx, r 'xi an d t'x ■ Assuming a large source-drain volt- 
age, such that all channels are equally populated, the 
correlation functions read: 

C21 = jjyTr {Ra — Ta} ; 

C 32 = ^Tr {i?^ (R B - T B ) + T\ (R' b -r B )}; 
C31 = ^Tr [rf A (R B - T B ) - r] (R' B - T' B )) 
+^Tr [r A t A (fJr B -r' B h B ) 

+Ur\{r B t' B -t B r' B )), (19) 

with Ra = t a ta, Ta = t A tA, etc. The second term in 
the expression for C31 arises from interference between 
the paths. In the single channel case, these results reduce 
to those of Eq. ([15]). 

An important observation can be made about the 
multi-channel case by considering that the scattering ma- 
trices preserve the channel-index, i. e. we essentially have 
M independent interferometers. In this case, the LG pa- 
rameter reads K = ]^Em=i^ m 'i where is the 



LG parameter for channel m. If we could tune by hand 
all the parameters of the scattering matrices, then the 
maximum violation of K max = | can be reached. How- 
ever, in an experiment, there will typically only be a 
few controllable parameters and this could make viola- 
tions hard to observe. Let us assume that we can adjust 
the parameters such that one of the K^> is maximised, 
m = 1, say. Whether we see a violation or not very much 
depends on what happens with the parameters of the 
other channels. If these parameters are all roughly sim- 
ilar to those of channel 1, then violations should still be 
observed. Generically, however, this will not be the case, 
and the for the other channels will take unrelated 

values in the range from —3 to |. The negative values 
are particulary troublesome as they will tend to over- 
whelm any positive contribution to the violation from 
other channels. This lack of controllability means that 
multi-channel geometries are best avoided if violations of 
the LGI are sought. 



III. DOUBLE QUANTUM DOT 
INTERFEROMETER 

The above MZI scheme functions by having the de- 
tectors remove electrons from the interferometer arms. 
In this section we study a second MZI realisation which 
leaves the electrons within the system and the effects of 
measurement are only felt through dephasing. This sec- 
ond set-up is shown in Fig. |4j As in the foregoing, the 
basic structure is of two (single-channel) leads that are 
joined at two beamsplitters. Beamsplitters between non- 
edge-channel leads can be realised by tunnel junctions, as 
in the recent experiments by Yamamoto et al^. In each 
arm of the interferometer there is a QD and alongside 
each QD is a QPC charge detector. When connected to 
a voltage supply, the current flowing through the QPCs 
serves as read-out of the occupation of their respective 
QDs. Note that although similar detectors we used in 
e.g. Refs. [6^ and [7(1 the way in which they are used here 
is different. 



A. Model 

We first consider the system without detectors. Our 
MZI model is related to that of, e.g., Refs. [n] and [72l. 
but with different leads. Far from the junctions, we 
describe the four leads as non-interacting Fermi reser- 
voirs with Hamiltonian H lcs = ^2 ^kaqc\ aq Ckaq with k 
the wavenumber of the electron, and where a = 1, 3 and 
q = ± specify the lead (we set h = 1 and ignore spin). 
We assume that there is but a single orbital of relevance 
in each dot and that the DQD system is in the strong 
Coulomb blockade regime, such that it is restricted to 
just three states: 'empty', |0); or with one excess elec- 
tron in either the upper or lower dots, |+) = gC_ |0) and 
|— ) = d_|0), respectively. Assuming the dot levels are 
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FIG. 4. (Color online) Sketch of a Mach-Zehnder interferom- 
eter with a quantum dot in each arm. The charge state of each 
QD can be monitored by the currents flowing through QPCs 
next to the dots. Here only the QPC monitoring the dot + 
is active, such that the correlation function to be measured is 
as Fig. [lb. 



detuned by an energy e from one another, the dot Hamil- 
tonian reads H$ = ^ e J2 q ld q d q . In the following, we set 
e = for simplicity. 

We assume that the effect of the beam-splitters is to 
modify the amplitudes with which the leads couple to 
the QDs. So, for example, an electron in lead 1+ tunnels 
into a superposition of upper and lower dot states, with 
the details of the superposition being determined by the 
scattering matrix s^. The tunnel Hamiltonian connect- 
ing lead and dots therefore reads 

H ? = E ( C Ii • S A ■ d + CL ■ sb ■ d + H.c.) , (20) 

k 

where sa,b are scattering matrices, assumed to be energy 
independent, d = (<f + ,<f_) is a vector of dot operators, 
and C\ a = {T a+ c\ a+ ,T a ^c\ a _) arc vectors of lead op- 
erators with tunnel matrix elements T as , also assumed 
to be energy- independent. The corresponding sequential 
tunnel rates are T aq = 2-n:\T aq \' 2 g aq , where g aq is the 
density-of-states of reservoir aq, also assumed constant. 

In the infinite-bias limit, the system can be described 
by a quantum master equation of Lindblad-formSi^. 
Let us introduce the super-operator notations Sf[d]p = 
dpS and A[d]p = — i {d^d, p^r— , and introduce the 
operators 

di q = \fi\~ q e q • sa ■ d; d 3q = \/T^e q ■ s B ■ d, (21) 

with unit vectors e + = (1,0) and e = (0,1). With 
introduction of counting fields Xa q to facilitate the cal- 
culation of current statistics (see e.g. Refs. l74l l79|). the 
X-resolved master equation for the DQD system reads 

p( X ) = -i [H S , P ]+Y. (e lXaq J[d ag ] - A[d aq }) p. (22) 

aq 

The QDs are monitored by QPCs in a set-up similar to 
the single dot in an interferometer in the experiment of 



Ref. HH. In including the detectors in our theory, we fol- 
low Gurvit a 58 ' 61 . When dot q is unoccupied, the Hamil- 
tonian for QPC Dq reads 

H Dq = J2 U3 k S qalsq a ksq + E a kLq a kRq + H.C, (23) 
ks k 

where oj® sq is the energy of an electron in state k on 
side s = L, R of the the QPC, and Vl q is the coupling 
amplitude between the two sides, (assumed energy inde- 
pendent). When dot q is occupied, we assume that this 
Hamiltonian is modified such that the coupling constants 
shift to different values, VL q — > fl q . In the limit of large 
bias across the QPC, the detector at location q gives rise 
to an extra Liouvillian 

W Dq (XD q ) = e l ^J[d Dq ] - A[d Dq ], (24) 

which adds to the DQD Liouvillian. Here, do q = 

V^l <?)(<? I + \Fh (! _ witn Iq the rate of elec " 

tron transfer through the QPC q when its dot is empty, 

and 7^ the rate when the dot is occupied. The count- 
ing field XD q here allows us to calculate the statistics 
of the detector currents. Microscopically, the rates are 
7 9 = 2ir\n q \ 2 p Lq p Rq V Dq and i q = 2ir\Vt' q \ 2 p Lq p Rq V Dq , 
with p sq the density of states of the QPC reservoir sq 
and Vu q the applied voltage. Detectors may be decou- 
pled or coupled from the MZI-QD system by adjusting 
the QPC voltages such that the differences between the 
amplitudes il q and Q' q is either zero (decoupled) or finite 
(coupled). Here, we only couple at most one detector to 
the system at a given time. Furthermore, we assume bal- 
anced detectors such that with the D+ detector coupled 
we have 7+ = 7, 7^ = 7', and 7_ = 7L = 0, and when 
the D— detector is coupled we have 7_ = 7, 7L =7', 
and 7+ = 7+ = 0, 



B. Current, correlation functions and probabilities 

Our approach to measuring the LGI with this set-up is 
similar to that with the quantum Hall edge-channels with 
the main exception being how C32 is obtained. We inject 
electrons into the '+'-channel of lead 1 and close the '1— ' 
channel: Ti + — > Tl and ri_ — > 0. For simplicity, we set 
the output rates equal: = T3- = Tr. 

To obtain C31, we switch off the QPC detectors and 
measure the output currents at 3±. Arranging the el- 
ements of the density matrix into a vector in the basis 

(00, ++, , + — , — h), the stationary state of the DQD 

system reads 

Artat = 5= (2r i? ,r L (l + cos6»A),rL(l-cos6'A), 

_^a/2 Tl sin A) _ e - 4 A /2p L sin6 / A ^ (25) 

with total width r = + T R . Here, we have assumed 
the same scattering matrices as in Eq. ([7]). The total 
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current flowing is (I) tot = (-0i+ = FlTr/T, which is 
divided between the output ports as 



(I) 



3± 



( 1 ± cos 8a cos Ob T cos 4> sin 8a sin ( 



Constructing C31 as in Eq. (|T2"j) we obtain 

C31 = cos6U cos 9b — cos (/) sin 9a sm9s, 



(26) 



which agrees with that of Eq. © . 

Next we can obtain C21 by turning on the QPC de- 
tectors one at a time. As shown in Ref. [Hg, the mean 
current flowing through the QPC can be used to extract 
the mean current flowing through the corresponding dot. 
With a detector coupled to dot q, the current through 
the detector is 



(I Dq ) 



(I)t 



2T, 



7 1 + 2- 



+ 7' + 9(7' - 7) cos 6a 



The current flowing through the QPC when the DQD is 
empty is (I% q ) = 7, such that the difference is 



(A/d 9 ) 



Dqj 



(I° Dq ) 



1±— !l(I) tot (I + q cos 9 A ), (27) 

2l R 



which is proportional to the probability that an electron 
takes the path 2q. Assuming balanced detectors, we ob- 
tain 

C 21 = { 4H-^4^os9 A , (28) 



(AIi 



(Air 



as in Eq. (jT2j) . 

Whereas these two correlation functions can be deter- 
mined with just mean-current measurements, to deter- 
mine C32 we need to consider current cross-correlations. 
Let us first imagine that we can measure the current 
through dot q. Then, in the limit Tl —> 0, such that 
there is only ever at most one electron in the interferom- 
eter at a given time, the zero-frequency noise correlator 



S3q'2q = 5/ ^{hq 1 , hq}) c 



where (. . .) c denotes the cumulant average, is propor- 
tional to the joint probability, P3 q '2q, that the electron 
travels through dot q and ends up at output 3q'. This 
result follows in the same way as in Ref. [l5|; the differ- 
ence here being that we correlate the position of a single 
electron in subsequent regions, as opposed to the corre- 
lation of two spatially-separate electrons. Measuring all 
four such correlators, we obtain the probabilities 



S?,„> 



3q'2q 



3q'2q 



r i r >->3r'2r 



(30) 



From these directly-obtained probabilities, we construct 
the idcal-ncgative-mcasuremcnt ones as 



pINM 
r 3q'2q 



P. 



3q' 



P3q'2q, 



where q = —q and the total probability at output 3q' is 
obtained from the currents 



L 3q) 



(32) 



These relations follow from charge conservation and the 
unidirectional nature of the transport. 

The QPC detectors couple not the current flowing 
though the dot, but rather to their occupations. In terms 
of the zero-frequency correlation function between cur- 
rent fluctuations in the QPC and those in one of the 3± 
ports, 
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3q'Dq 



i J di({I 3q ,J Dq }) ( 



the required probabilities read 



P. 



3q'2q 



(Al£> q )S3 q 'p q 



(33) 



(34) 



This can be understood as follows. Whereas Ss q '2 q corre- 
lates two delta-function peaks corresponding to the pas- 
sage of the electron through the regions 2 and 3, S^q'Dq 
correlates a delta-function in region 3 with a signal of fi- 
nite duration in region 2, which corresponds to the finite 
time for which the dot is occupied. This mean occupa- 
tion time is proportional to the inverse of mean current 
through the dot, which can be obtained (up to a pro- 
portionality constant) from the mean detector current 

(AI D q). 

Calculating these probabilities, we find that in the 
limit — ^ 0, the third correlation function reads 



C32 = cos6» B , (35) 

in accordance with Eq. (|10p . Since, in the — > limit, 
all three correlation function are identical with their ideal 
counterparts, the LCI for this set-up is identical to that of 
Eq. (fTTj) . In the way that we have described the QPC de- 
tectors here, it docs not make any difference whether we 
calculate K using the ideal negative measurement prob- 
abilities or the direct ones since, in our theoretical de- 
scription, the QPC detectors act as ideal detectors and 
(29) only influence the system through their dephasing effect. 
Experimentally, the ideal negative measurement protocol 
should be used, and actually, the comparison between the 
case with ideal negative measurement and that without 
would give an interesting method for studying to what 
extent the QPC measurements are non-invasive. Let us 
just add that, whilst the above results were derived in the 
symmetric case with T^/Tr- = 1 and with balanced 
detector rates, if these ratios are unequal but known, 
then the difference can be accounted for by weighting 
the terms in the correlation functions accordingly. 



C. Dephasing 



A simple way to include the effects of dephasing in 
(31) this model is to "leave the detectors switched on" when 



8 



calculating C31. With empty and occupied rates 7dephase 
and 7d phasc the measured K function has the form of 
Eq. (fT6|) . with the function /(A) replaced by 



(V V-/ 2 

\ 'dephase/ 



("Tdephase)) 



1/2 



2r, 



(36) 



/ = 1 is the ideal no-dephasing case, and / = gives the 
classical limit. To obtain strong violations of the LGI, 
therefore requires that the difference in rates 7d op hasc an< ^ 
7dephase is small compared with the tunnel rate T R . 



D. Detection errors 

Just as the direct relation between the Bell inequal- 
ities and noise measurements of, e.g., Refs. [l5| and [H 



relies on the weak-tunnel limitiir— , so it is here that our 
measurements are only isomorphic with those required 
by the LGI in the F^ — >• limit. Away from this limit, 
there exists the possibility that our measurements mis- 
takenly correlate subsequent electrons, rather than the 
same electron with itself. 



The LGI quantity can be calculated using the currents 
and zero-frequency noise, as described above, away from 
the — > limit to assess the error. Assuming for sim- 
plicity that the detector is faster than the system dynam- 
ics 7' — 7 3> Vl/r (although the general case can easily 
be investigated too), we obtain for the LG correlator 



K' 



{ [- (T L + T R f + T L (T L + 3T R ) cos 2 6 A ] cos 8 B 



(L L - Y R )Y R 

+ -r fl )r fl [2cos6» A sin 2 (6» B /2) +cos0sin6» A sin0 B ]} 



(37) 



This expression is again maximised with cos<^> = 1, but, 
unlike the L^ — > case, the maximizing angles 8 a and 
6b are not equal. If we assume that Tl/Y r <C 1, we can 
expand to leading order (7' — 7 — > 00) to obtain 

K" = K + 3^- cosOb sin 2 9 A + O^^y^j , 

where K is the L^ — >• value. We can also calculate 
the corresponding quantity in the classical limit (this we 
do by calculating C31 in limit (7dc P hasc ~ 7dc P hase) -» 
00). In this case, we obtain C$\ = cos 8 a cos 8b, and the 
expansion of B' = C21 + C32 — C%\ for small L^ gives 

B' = B + 3^ cosdB sin 2 8 A + 0^^y^j , (38) 

where B = cos 9 a + cos 9b — cos9acos9b is the classi- 
cal value in the ideal case which, when maximised gives 
Bmux = 1, the bound of Eq. ([T]). Maximising B' over 
the angles, we obtain a value bigger that unity. To low- 
est order then, classical and quantum LG correlators are 
affected in the same way. Fig. [5] shows the maximum 
values of both quantum and classical correlators. 

Thus, assuming that we know the ratio of T^/T R 
from current and noise measurements, the effects of a 
finite tunneling rate F^, can be included in assessment 
of whether LGI is violated or not. The conservative ap- 
proach is say that the quantity (K' mllx — |) represents a 
systematic error in the measurement, and assuming that 
this error works against us, we can only conclude that we 
violate the LGI when the measured value of K exceeds 



unity by an amount equal to this error. Alternatively, one 
can say that since one knows how the classical bound be- 
haves at finite T/T R , we can simply use B' max of Eq. (|38[) 
as a bound. However, providing that we are in the cor- 
rect operating limit of Tl/T r <C 1, these modifications 
will be very small, such that whether they are taken into 
account or not will only effect the question of violation 
in marginal cases. 



IV. CONCLUSIONS 

We have considered the violation of the LGI in MZ 
inteferometer geometries. The key to the violation is a 
combination of the interference at the second beamsplit- 
ter and the inhibition of this interference by the measure- 
ment process. In the two proposals we have considered 
this inhibition occurs in two different ways. In the first 
realisation, we physically interrupt transmission through 
one of the arms of the MZI, obviously preventing inter- 
ference. On the other hand, in the DQD proposal, the 
detectors act in a more traditional way and introduce 
dephasing between the paths. 

In this MZ geometry both the state of the electron 
and measurement time are mapped onto real-space co- 
ordinates — the qubit states Q = ± are physically sep- 
arate paths, and the regions within the interferometer 
correspond to different time instances. This mapping 
has several advantages for seeking a violation of the LGI 
in transport. The mapping of the time-coordinate means 
that we do not need to make time-resolved correlation 
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I U-^ i I i I i I i I 

0.05 0.1 0.15 0.2 

r /r 

l' r 

FIG. 5. (Color online) The maximum value K' max (blue solid 
line) and the corresponding classical value B^ax (black dashed 
line) away from the Tl — > limit. Both are higher than their 
ideal Tl — > values. A symmetric system was assumed and 
the fast detector limit (7' — 7) /Tr — > 00 taken. 



measurements. All the measurements required here are 
either mean stationary currents or zero-frequency noise 
correlators. Furthermore, the spatial separation of the 
qubit degrees of freedom facilitates the realisation of ideal 
negative measurement, since it is relatively easy to cou- 
ple to just one of the qubit states when they are spatially 
distinct. In this respect, increasing the separation of the 
detector arms should decrease the plausibility of claims 
that detection in one arm is, from a macro-realist point- 
of-view, influencing the other. 

The general principles described here can easily be ex- 
tended to further systems. Within transport, for exam- 
ple, our second scheme could also be realised with an 



edge-channel MZI plus QPC detector channel without 
the quantum dots^ r 49 ' 51 . An alternative setting for the 
realisation of our first scheme might be the flying qubit 
experiment of Ref. H3, which is essentially a MZI away 
from the quantum Hall regime. Two challenges are ob- 
vious with this realisation. Firstly, the leads reported in 
the experiment have multiple channels, which potentially 
gives rise to the problems discussed in section Hi CI The 
second problem is that of backscattcring at the beam- 
splitters and detectors, which has (justifiably) been ne- 
glected here but probably can not be eliminated in set- 
ups such as that of Ref. [57|. 

Applications away from electronic transport are also 
possible. The application of the first scheme in optics is 
obvious but the notion of the qubit state is predicated 
on the source being a single-photon source. So whilst a 
classical wave might also exceed the right-hand-side of 
Eq. ([lj) , this would not constitute a violation of the LGI, 
as it represents an application of the concepts outside 
their proper realm of definition (i.e., non-dichotomic ob- 
servables). Going further, the same principles could be 
used to test the LGI with electrons in free space, neu- 
trons, atoms and molecules, all of which have had in- 
terference experiments in the MZI geometry conducted 
on the m 80 ' 81 . Of these, molecules offer the most excit- 
ing prospect, as there the nature of the coherence being 
tested could potentially be macroscopic, in line with the 
original goals of Ref. y. 
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